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Abstract. The inductive dimension dim(G) of a finite undirected graph G 
is a rational number defined inductively as 1 plus the arithmetic mean of 
the dimensions of the unit spheres dim(S'(a;)) at vertices x primed by the 
requirement that the empty graph has dimension —1. We look at the dis- 
tribution of the random variable dim on the Erdos-Renyi probability space 
G(n,p), where each of the n(n — l)/2 edges appears independently with prob- 
ability < p < 1. We show that the average dimension dn{p) = Ep,n[dim] 
is a computable polynomial of degree n{n — l)/2 in p. The explicit formu- 
las allow experimentally to explore limiting laws for the dimension of large 
graphs. In this context of random graph geometry, we mention explicit for- 
mulas for the expectation Ep,n[x] of the Euler characteristic x, considered 
as a random variable on G(n,p). We look experimentally at the statistics of 
curvature K(v) and local dimension dim(^') = 1 + dim(S'(i;)) which satisfy 
x{G) = ^y^v dim(G) = ^ J^^^y dim('u). We also look at the 

signature functions f(p) = Kp[dim\, g(p) = Ep[x] and matrix values functions 
Av,w(p) = Covp[diin(v),diin(w)], By, w{p) = Cov[K (v) , K (w)] on the proba- 
bility space G(p) of all subgraphs of a host graph G = (V, E) with the same 
vertex set V, where each edge is turned on with probability p. 

1. Dimension and Euler characteristic of graphs 

The inductive dimension for graphs G = (V^E) is formahy close to the Menger- 
Uhryson dimension in topology. It was in [7 defined as 

dim(0) = -l,dim(G) = 1 + XI , 

where S{v) = {w e V \ {w,v) e E },{e = {a,b) e E \ {v,a) e E,{v,b) e E } 
denotes the unit sphere of a vertex v e V. The inductive dimension is use- 
ful when studying the geometry of graphs. We can look at the local dimen- 
sion dim('L') = 1 + dim(S'(i')) of a vertex which is a local property like "degree" 
deg{v) = OTd{S{v))^ or "curvature" K{v) defined below. Dimension is a rational 
number defined for every finite graph; however it is in general not an integer. 

The dimension is zero for graphs Pn of size 0, graphs which are completely dis- 
connected. It is equal to n — 1 for complete graphs Kn of order n, where the size 
is (2) = n(n — l)/2. Platonic solids can have dimension 1 like for the cube and 
dodecahedron, it can be 2 like for the octahedron and icosahedron or be equal to 
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Figure 1 . This figure shows a particular graph G with 20 vertices 
and Euler characteristic x{G) = 1. The left side shows the graph 
with curvatures 

{0, -1, -2, 6, 3, -5, -2, -4, 3, 0, -4, -1, 0, 2, 1, 1, 1, 2, 3, 3}/6 

which sum up to the Euler characteristic 1 by Gauss-Bonnet. The 
right figure shows the same graph with local dimensions 

{2, 5/3, 2, 0, 1, 8/5, 7/4, 3/2, 1, 1, 11/6, 5/3, 1, 2, 2, 2, 2, 2, 1, 1 } 

which average to the dimension of the graph dim(G) = 1801/1200. 



3 like for the tetrahedron. The 600 cell with 120 vertices is an example of a three 
dimensional graph, where each unit sphere is a two dimensional icosahedron. Fig- 
ure [3] illustrates two Archimedean solids for which fractional dimensions occur in 
familiar situations. All Platonic, Archimedean and Catalan solids are graph the- 
oretical polyhedral a finite truncation or kising process produces two-dimensional 
graphs. 

Computing the dimension of the unit spheres requires to find all unit spheres 
for vertices in the unit sphere S{x) of a vertex and so on. The local dimension 
dim(x) = 1 + dim(5'(x)) satisfies by definition 

(1) dim(G) = E ^i^(^) • 

There is an upper bound of the dimension in terms of the average degree deg(G) = 
JV] deg('y) of a graph: dim(G) < deg(G) — 1 with equality for complete graphs 
Kn. The case of trees shows that dim(G) = 1 is possible for arbitrary large |V| or 
deg(G). 

An other natural quantity for graphs is the Euler characteristic 



(2) 



00 

X(G) = ^(-l)^fe , 

k=0 
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where Vk is the number of Kk-\-i subgraphs in G. We noted in ^ that it can be 
expressed as the sum over all curvature 

k=0 

where Vk{x) is the number of Kk-\-i subgraphs in the sphere S{x) at a vertex x. As 
in the continuum, the Gauss-Bonnet formula 

(3) x{G)=^K{x) 

xev 

relates a local quantity curvature with the global topological invariant x- For 
example, for a graph without triangles and especially for trees, the curvature is 
K{v) = 1 — deg{v)/2. For graphs without 4 cliques and especially two dimensional 
graphs K{v) = 1 — deg{v)/2 — size{S{v)) /3. For geometric graphs, where each 
sphere S{v) is a cyclic graph, K{v) = 1 — deg{v)/6. For the standard Petersen 
graph P5,2, the dimension is 1, the local dimension constant 1, the Euler character- 
istic —5 and the curvature is constant —1/2 at every vertex. The Petersen graph 
Pg^s has dimension 4/3 and Euler characteristic —6. There are 9 vertices with cur- 
vature — 1/2 and 9 with curvature —1/6. The sum of curvatures is —6. 

The Gauss-Bonnet relation (|3| is already useful for computing the Euler char- 
acteristic. For inhomogeneous large graphs especially, the Gauss-Bonnet-Chern 
formula simplifies in an elegant way the search for cliques in large graphs. An other 
application is the study of higher dimensional polytopes. It follows immediately 
for example that there is no 4-dimensional polytope - they are usually realized as a 
convex set in - for which the graph theoretical unit sphere is a three dimensional 
600 cell: the Euler characteristic is 2 in that dimension and the curvature would 
by regularity have to be 2/|V|. But curvature of such a graph would be constant 
and since the 600 cell has 120 vertices, 720 edges, 1200 faces and 600 chambers, the 
curvature were constant 

V_i Vo Vi V2 Vs 1 120 720 1200 600 , 

^ = — -Y + y-T + T = i — + — + ^ = 1 

which would force |y| = 2 and obviously does not work. While such a result could 
certainly also be derived also with tools developed by geometers like Schlafii or 
Coxeter, the just graph theoretical argument is more beautiful. It especially does 
not need any ambient space realization of the polytope. 

We have studied the Gauss-Bonnet theme in a geometric setting for d-dimensional 
graphs for which unit spheres of a graph satisfy properties familiar to unit spheres 
in R^. In that case, the results look more similar to differential geometry [8]. 

The curvature for three-dimensional graphs for example is zero everywhere and 
positive sectional curvature everywhere leads to definite bounds on the diameter of 
the graph. Also for higher dimensional graphs, similar than Bonnet-Schoenberg- 
Myers bounds assure in the continuum, positive curvature forces the graph to be of 
small diameter, allowing to compute the Euler characteristic in finitely many cases 
allowing in principle to answer Hopf type questions about the Euler characteristic 
of finite graphs with positive sectional curvature by checking finitely many cases. 
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Figure 2. The Petersen graph P(12,4) of dimension 4/3 and 
Euler characteristic —8 can be seen to the left. The curvatures 
take 12 times the value —1/2 and 12 times the value —1/6. To the 
right, we see the Turan graph T(13,4) of dimension 3 and Euler 
characteristic —23. The curvatures take values —2 nine times and 
—5/4 four times. 




Figure 3. The dimension of the small rhombicosidodeca- 

hedron is 3/2: each point has dimension 3/2 because each unit 
sphere is a graph with 4 vertices, where two points have dimension 
1 and two points have dimension 0. The unit sphere has dimen- 
sion 1/2. The dimension of the truncated dodecahedron is 5/3 
because every unit sphere has dimension 2/3. 

Such questions are often open in classical differential geometry but become a finite 
combinatorial problem in graph theory. Obviously and naturally, many question in 
differential geometry, whether variational, spectral or topological can be asked in 
pure graph theory without imposing any additional structure on the graph. Dimen- 
sion, curvature and Euler characteristic definitely play an important role in such 
quests. 

What are the connections between dimension and Euler characteristic and cur- 
vature? We don't know much yet, but there are indications of more connections: 
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all flat connected graphs we have seen for example are geometric with uniform con- 
stant dimension like cyclic graphs or toral graphs. An interesting question is to 
study graphs where Euler characteristic is extremal, where curvature (also Ricci 
or scalar analogues of curvature) is extremal or where curvature is constant. So 
far, the only connected graphs with constant curvature we know of are complete 
graphs Kn with curvature 1/n, cyclic graphs Cn with curvature 0, discrete graphs 
Pn with curvature is 1, the octahedron O with curvature 1/6, the icosahedron / 
with curvature 1/12, higher dimensional cross poly topes with 2n vertices and cur- 
vature l/(2n), the 600 cell 120 vertices, where each unit sphere is an icosahedron 
and which is "flat" 

1 2 3 4 1 2 3 4 

like for any three dimensional geometric graph [8 , twisted tori i^n,m with curvature 
as well as higher dimensional regular tesselations of tori. We could not yet con- 
struct a graph with constant negative curvature even so they most likely do exist. 
We start to believe that geometric graphs - for which local dimension is constant 
like the ones just mentioned - are the only connected constant curvature graphs. 

We look in this article at connections with random graph theory (TJ |10], an 
area of mathematics which has become useful for the study of complex networks 
[H m [HI E]- The emergence of interest in web graphs, social networks, neural 
networks, complex proteins or nano technology makes it an active area of research. 



2. Random subgraphs of the complete graph 

We inquire in this section about the dimension of a typical graph in the prob- 
ability space G(n, 1/2), where each edge is turned on with probability 1/2. To 
investigate this, we can look at all possible graphs on a fixed vertex set of cardi- 
nality n and find the dimension expectation by computing the dimension of each 
graph and adding this up. When we counted dimensions for small n by brute 
force, we noticed to our surprise that the sum of all dimensions of subgraphs is an 
integer. Our limit for brute force summation was n = 7, where we have already 
2^^ = 2^097^152 graphs. For the next entry n = 8, we would have had to check 128 
times more graphs. 

Note that as usual, we do not sum over all subgraphs of but all subgraphs of 
Kn for which \V\ = vq = n. While this makes no difference for dimension because 
isolated points have dimension 0, it will matter for Euler characteristic later on, 
because isolated points have Euler characteristic 1. For for example, there are 3 
graphs with dimension 1, one graph with dimension 2 and 3 graphs with dimension 
2/3. The sum of all dimensions is 7. For i^4, the sum of dimensions over all 2^ = 64 
subgraphs is 75: there are 22 subgraphs of dimension 1, there are 12 of dimensions 
5/3 and 3/4 each, 6 of dimensions 1/2 and 2 each, 4 of dimension 3/2 and one of 
dimension and 3 each. Also for K5, dimension 1 appears most with 237 subgraphs 
followed with 120 graphs of dimension 22/15. For Kq already, we have two integer 
champions: dimension 1 appears for 3^981 subgraphs and dimension 2 for 2^692 
subgraphs. The total sum of dimension is 53^275. 
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Theorem 2.1 (Average dimension on G(n, 1/2)). The average dimension on G{n, 1/2) 
satisfies the recursion 

dn+1 = 1 + X] ' 

where = —1 is the seed for the empty graph. The sum over all dimensions of all 
order n subgraphs of Kn is an integer. 

Proof. Let g{n) = 2(3) be the number of graphs on the vertex set V = {1, n }, 
and let f{n) the sum of the dimensions of ah subgraphs of the complete graph 
with n vertices. We can find a recursion for f{n) by adding a (n + l)'th vertex 
point X and then count the sum of the dimensions over all subgraphs of iCn+i 
vertices by partitioning this set of subgraphs up into the set Yj^ which have k edges 
connecting the old graph to the new vertex. There are (^) possibilities to build 
such connections. In each of these cases, the unit sphere S{x) is a complete graph 
of k vertices. We get so the sum Xl^=o(/fc ^ dimensions of such graphs 

because we add 1 to each of the (^) cases. From this formula, we can see that the 
sum of dimensions is an integer. The dimension itself satisfies the recursion 



^/,+i = 2(s)-m^(4 + i) 



k=0 



With (^) - (^+^) = -n and using 2"^ ^"^^^ il) = we get the formula. □ 
Remarks. 

1. We see that (in+i is 1 plus a Bernoulli average of the sequence {d^ 1^=0- 



2. Theorem 2.1 will be generalized in Theorem 4.1 from p = 1/2 to general p. 



3. Average dimension on classes of graphs 

We can look at the average dimension on subclasses of graphs. For example, 
what is the average dimension on the set of all 2-dimensional graph with n ver- 
tices? This is not so easy to determine because we can not enumerate easily all 
two-dimensional graphs of order n. Only up to discrete homotopy operations, the 
situation for two dimensional finite graphs is the same as for two dimensional man- 
ifolds in that Euler characteristic and the orientation determines the equivalence 
class. 



We first looked therefore at the one-dimensional case. Also here, summing up 
the dimensions of all subgraphs by brute force showed that the sum is always an 
integer and that the dimension is constant equal to 3/4. This is true for a general 
one dimensional graph without boundary, graphs which are the disjoint union of 
circular graphs. 

Theorem 3.1 (3/4 theorem). For a one- dimensional graph G = {V^E) without 
boundary, the average dimension of all subgraphs of G is always equal to 3/4. The 
sum of the dimensions of subgraphs H of G with the same vertex set V is an integer. 

Proof. A one-dimensional graph without boundary is a finite union of cyclic graphs. 
For two disjoint graphs Gi,G2 with union G = Gi U G2, we have dim(G) = 
(dim(Gi)|Gi| + dim(G2)|G2|)/|G|. It is therefore enough to prove the statement 
for a connected circular graph of size n > 4. This can be done by induction. We 
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can add a vertex in the middle of one of the edges to get from Cn to Cn+i- The 
dimensions of the other vertices do not change. The new point has dimension 
with probability 1/4 and 1 with probability 3/4. Since the smallest one dimensional 

graph has 4 nodes, and is already a multiple of 4, the sum dimensions of all 
subgraphs is an integers. □ 



Remarks. 

1. We will generalize this result below and show that for a one dimensional circular 
graph, the expected dimension of a subgraph is p{2—p). The 3/4 result is the special 
case when p = 1/2. The result is of course different for the triangle Ks = C3, which 
is two dimensional and for which the expected dimension is p{2 — p ^ p^). We will 
call the function f{p) = p{2 —p) the signature function. It is in this particular case 
the same for all one dimensional graphs without boundary. 

2. Is the sum of dimensions of subgraphs of a graph with integer dimension and 
constant degree an integer? No. Already for an octahedron, a graph G of dimension 
2 which has 2^^ = 4096 different subgraphs, a brute force computation shows that 
the sum of dimensions is a = 15424/3 and that the average dimension of a subgraph 
of G is a/2^^ = 1.25521. The unit ball Bi(v) in the octahedron is the wheel graph 
W4 with four spikes in which the sum of the dimensions is 284 and the average 
dimension is 284/256 = 1.10938. 

3. We would like to know the sum of all dimensions of subgraphs for flat tori, finite 
graphs for which each unit disc is a wheel graph Wq. Such tori are determined by 
the lengths M, K of the smallest homotopically nontrivial one dimensional cycles 
as well as a "Dehn" twist parameter for the identification. Already in the smallest 

/48\ 

case M = = 4, the graph has 48 edges and summing up over all 2^2^ = 3.6 • 
10^^^ possible subgraphs is impossible. It becomes a percolation problem [5 . It 
would be interesting to know for example whether the dimension signature functions 
/k,m(p) = Ep[dim] have a continuous limit for min(i^, M) 00. 

4. For the wheel graph Wq^ the unit disc in the fiat torus, the average dimension 
is ( 159368/35) /2^^ = 1.1116... For a flat torus like Ts^s we can not get the average 
dimension exactly but measure it to be about 1.3. We observe in general however 
that the dimension depends smoothly on p. 

5. Brute force computations are not in vain because they allow also to look at the 
distribution of the dimensions of graphs with n vertices. Since we do not yet have 
recursive formulas for the higher moments, it is not clear how this behaves in the 
limit n ^ 00. 

The first sums of dimensions are Si = 0, 52 = 1, 53 = 7, 54 = 75, 55 = 1451, sq = 
53275, 57 = 3791451, sg = 528662939, sg = 145314459035, sio = 79040486514843, 
and 5ii = 85289166797880475. On a graph with 5 vertices for example, the sum of 
all dimensions over all the 2^^ subgraphs is 1451. We could not possibly compute 
5ii by brute force, since there are just too many graphs. 

Favoring integer dimensions for concrete or random graphs is a resonance phe- 
nomena of number theoretical nature. Whether it is a case for "Guy's law of small 
numbers" disappearing in the limit n ^ 00 remains to be seen. 
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Figure 4. The distribution of the random variable dim on the 
finite probabihty space G(7, 1/2) of ah graphs with vertex cardinal- 
ity 7. The left picture shows the cumulative distribution function, 
the right the density function. Graphs with integer dimension ap- 
pear to have higher probability than graphs with fractional dimen- 
sion. This is a phenomenon, we also see in concrete networks, like 
social graphs or computer networks. The function dim takes 245 
different dimension values on G(J^ 1/2). The dimension 2 appears 
for 4^146^583 graphs, the next frequent dimension is 1 appears in 
99^900 graphs. Only in third rank is dimension 55/42 is a fraction. 
It appears for 55^440 graphs. 



I I I I ll lllll 1 4 I I llllll J . , 1 . 

Figure 5. The dimension statistics on G(5,l/2) and G(6, 1/2) 
already indicates a preference for integer dimension of subgraphs. 

4. The dimension of a random p-percolating graph 

We generalize now the probability measure on the space of all graphs with n 
elements and switch each edge on with probability < p < 1. This is the classical 
Erdos-Renyi model [3]. With the probability measure Pp on X^, the probability 
space is called G{n,p). For this percolation problem on the complete graph, the 
mean degree is E[deg] = np. 
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The following result is a generalization of Theorem 2.1 , in which we had p = 1/2. 
It computes = En, p [dim]. 

Theorem 4.1 (Average dimension on G{n,p)). The expected dimension Ep[dim] 
on G{n^p) satisfies 

where do = —1. Each dn is a polynomial in p of degree (2). 

Proof. The inductive derivation for p = 1/2 generalizes: add a n + I'th point P and 
partition the number of graphs into sets Y^, where P connects to a /c-dimensional 
graph within the old graph. The expected dimension of the new point is then 



n . 

fe=i ^ 



and this is also the expected dimension of the entire graph. This can be written as 
d{n + I) = j2 - p)"-'=rf(fc) + - 

k=l ^ ^ k=l ^ ^ 



which is 

n / 

d(n + 1) = ^ r - pr->^d{k) + 1 - (1 - p)"d(o) - (1 - py 



k=0 

which is equivalent to the statement. □ 

Again, if we think of the vector d = (do, • • • , <^n) as a random variable on the 
finite set {0, 1, . . . , n } then dn-\-i is 1 plus the expectation of this random variable 
with respect to the Bernoulli distribution on {0, 1, . . . , n } with parameters n and p. 



Lets look at the first few steps: we start with —1, where the expectation is 

— 1 and add 1 to get 0. Now we have (—1,0) and compute the expectation of this 
(— l)p^(l— p)^+Op^(l— p)^ = p—1. Adding 1 gives p so that we have the probabilities 
(— l,0,p). Now compute the expectation again with {—l)p^{l — p)^ + 2p^{l — p)^ ■ 
0+p^(l — Adding 1 gives the expected dimension of ds{p) = 2p — on a 
graph with 3 vertices. We have now the vector (—1, 0,p, 2p — p^ To compute 
the dimension expectation on a graph with 4 vertices, we compute the expectation 
(-1)/(1 -pf + 3p^(l - • + ?>p^(l -pf -p^ lp^(l -pf ■ (2p-p^ + p^) = 

— 1 + 3 • p — 3 • + 4 • — — + and add 1 to get the expected dimension 
d4^{p) = 3 ■ p — 3 - p'^ -\- 4: ■ p^ — p"^ — p^ -\- p^ on a graph of 4 elements. Here are the 
first polynomials: 

d2{p) =P 
ds{p) = 2p-p2 

d^ip) = 3p- 3p^ + Ap^ -p^ -p^ ^p^ 

d^lp) =4p- 6p'^ + lOp^ - 5p^ - 3p^ + 5/ - / - / + p^^ 

The expected number of cliques Kj^ of the complete graph Kn is {^)p^^^ • By ap- 
proximating the Binomial coefficients (^) with the Stirling formulas, Erdos-Renyi 
have shown (see Corollary 4 in [4 ), that for p < n~'^^^^~^\ there are no K]^ sub- 
graphs in the limit n ^ 00 which of course implies that then the dimension is 
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Figure 6. The expected dimension of a graph as a function of p 
is monotone and is shown here for n = 1, . . . , 20. We see the graphs 
of the functions dn{p)/{n — 1) for n = 1, 2, . . . 20. A division by 
(n — 1) produces a scahng so that dn{l) = 1. We experiment ahy 
see that they satisfy a power law: there are functions constants 
a = a{p) and c = c{p) for which with the scahng law dn{p) ~ c/n^ 
holds in the limit n ^ oo. Especially dn{p)/{n — 1) is monotone in 
n for every p. This needs to be explored more. 




Figure 7. Three random graphs in which the percolation prob- 
ability p parameter is chosen to have an expected dimension 1, or 
2 or 3. 

< {k — 1) almost surely. For p < n~^/^, there are no triangles in the limit and 
dim(G) < 1 for p < n~^/^, there are no tetrahedra in the limit and dim(G) < 2. 

Is there a threshold, so that for p < the expectation of dimension converges? 
We see experimentally for any p that E[dim] ~ c(jp)/n'^^^\ where c, a depend on p. 

Finally lets look at a generalization of the 3/4 theorem: 

Theorem 4.2 (Expected dimension of a subgraph of a one dimensional graph). 
The expected dimension Ep[dim] on all subgraphs of any one- dimensional graph 
without boundary is p{2 — p). 
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Figure 8. The sum Sn{k) of dimensions of size-k and order-n 
subgraphs of a complete graph is an integer too. The sequence 
Sn{k) appears to converge to a Gaussian distribution. This is not 
surprising, when considering the de Moivre-Laplace hmit theorem 
and that dimension is highly correlated with the size of the graph. 
We see in experiments, that this distribution limit is achieved also 
for subgraphs of other graphs like circular graphs. 

Proof. Proceed by induction. Add an other vertex in the middle of a given edge. 
The expectation of the dimension of the remaining points does not change. The 
expectation of the dimension of the new point is l-p^ + 2p(l— p)+0(l— p)^ = p{2—p) 
because the dimension of the point is one if both one of the two adjacent edges are 
present and if none is present. □ 

Remark. More interesting is to introduce the variable q = p = 1 and write the 
expectation as a sum 



where ak^n is the expectation of the dimension on the probability space G{p) of all 
graphs (V, £^/e), where Ej^ runs over all subsets of the vertices of the host graph G. 
For example, for the host graph Cs we have 

E[dim] = 2(7^p + 13(7 V + 36(7 V + 55(7 V + 50(7 V+ 27(7 V + 8(7 V+/ =p{2-p) 

from which we can deduce for example that 27/28 is the expected dimension if two 
links are missing in a circular graph Cs or 13/28 if two links are present and that 
2/8 is the expected dimension if only one link is present. The polynomial has the 



5. The Euler Characteristic of a random graph 

Besides the degree average deg and dimension dim, an other natural random 
variable on the probability space G{n^p) is the Euler characteristic x defined in 



For p = 1/2, we sum up the Euler characteristics over all order n sub graphs of 
the complete graph of order n and then average. The list of sums starts with 




form 



2p(7^-^ + (2n - 3)(7''" V + n^g""" V + • • • + nqp^'^ + p' 



Equation 



1,3,13,95,1201,25279,... 
leading to expected Euler characteristic values 

1, 1.5, 1.625, 1.48438, 1.17285, 0.771454, 0.34842, -0.0399132 . 
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Figure 9. The total number of triangles K^, in all subgraphs of 
Kn is 2(2) (^)/8. To the left, we see ah the 2^ = 64 subgraphs 
of i^4. There are 32 triangles in total so that £['^3] = 1/2. To 
the right we see the statistics of the Euler characteristic on the 
probability space G(7, 1/2). This is the last row of Table [1] The 
most frequent case is x = 1 followed by x = 0. 



Already for n = 5, there are some subgraphs with negative Euler characteristics. 
For n = 8, the expectation of x has become negative for the first time. It will do so 
again and again. These first numbers show by no means any trend: the expectation 
value of X will oscillate indefinitely between negative and positive regimes and grow 
in amplitude. This follows from the following explicit formula: 

Theorem 5.1 (Expectation of Euler characteristic). The expectation value of the 
Euler characteristic on G{n,p) is 

k=l ^ ^ 

Proof. We only need expectation of the random variable Vk on G{n^p). But this is 
well known: (see [1 ). We have 

EH 

This later formula is proven as follows: if Sk is the set of all k-cliques in the 
graph. Now count the number n{s) of times that a k-clique s appears. Then 
E['i;/e]2(2) = XlsG^fc ^('^) ~ because n{s) is constant on Sk- We especially 

have E[vo] = n^E[vi] = {^)p is the expected value of the number of edges, and 
E[v2] = (3)^^ is the expected number of triangles in the graph and E[vs] = (4)^^ 
is the expected number of tetrahedra K4 in the graph. □ 

Remarks. 

1. The Euler characteristic expectation E^^p[x] as a function of n oscillates between 

different signs for n ^ 00 because for each fixed /c, the function n ^ (2) /2^'^^ dom- 
inates in some range than is taken over by an other part of the sum. 

2. We could look at values of Pn for which the expectation value of the dimension 
gives E^^p[dim] = 2 and then look at the limit of the Euler characteristic. 

3. The clustering coefficient which is 3v2/a2 where a2 are the number of pairs of 
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1 





n=7 35 


420 6321 


54985 


307475 


734670 


839910 


133693 


17206 


2170 


245 


21 


1 



Table 1. The Euler characteristic statistics of subgraphs of Kn- 
We see that x = 1 is the most frequent situation in ah cases n = 1 
to n = 7. 



adjacent edges is a random variable studied already. It would be interesting to see 
the relation of clustering coefficient with dimension. 

4. The formula appears close to Ylk=ii~^)^ {T)P^ which simplifies to 1 — (1 — p)^ 
and which is monotone in p. But changing the p^ to p^'^) completely changes the 
function because it becomes a Taylor series in p which is sparse. We could take 
the sequence E^^p[x] 7(^/2) keep the functions bounded, but this converges to 
for j9 < 1. A natural question is whether for some we can achieve that 
En,p(n)[x] converges to a fixed prescribed Euler characteristic. Since for every n, 
we have E^^o[x] = ^ E^^o[x] = 1 because of continuity with respect to 
we definitely can find such sequences for prescribed X ^ 1- 

5. Despite correlation between the random variables Vk and diff'erent expectation 
values, its not impossible that the distribution of the random variable X^(G) = 

— E[x{G)])/a{x{G)) on G{n,p) could converge weakly in the limit n ^ oc 
if a{x) is the standard deviation. We would like therefore to find moments of x 
on G{n^p). This looks doable since we can compute moments and correlations 
of the random variables Vk- The later are dependent: the conditional expectation 
E[t'3 = l|t'2 = ] for example is zero and so diff'erent from E[t'3 = 1] •E[t'2 = ] > 0. 

6. When measuring correlations and variance, we have no analytic formulas yet 
and for n > 8 we had to look at Monte Carlo experiments. Tests in cases where 
we have analytic knowledge and can compare the analytic and experimental results 
indicate that for n < 15 and sample size lO'OOO the error is less than one percent. 



di{p) 


= 1 




d2{p) 


= 2 


-P 


dsip) 


= 3 


-Sp^p^ 


d4.{p) 


= 4 


- 6p + 4p^ - p^ 




= 5 


- 10p-\-10p^ - bp^ 



Table 2. Polynomials dn{p) which give the average Euler charac- 
teristic Ep[x] on the probability space G{n,p). 



6. Statistical signatures 

For any graph G = (V, E) called "host graph" we can look at the dimension and 
Euler characteristic signature functions 

f{p) = Ep[dim(if) ], g{p) = Ep[x(//) ] 
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Figure 10. The average Euler characteristic as a function of 
n oscihates. The intervals, where the behavior is monotone gets 
larger and larger however, the amplitudes grow too. Also the p 
dependence is interesting. The right picture shows the average 
Euler characteristic as a function of p for a fixed order n = 52. 
The function is n — 1 at p = and 1 at p = 1. 

which give the expected dimension and Euler characteristic on the probability space 
G{p) of all subgraphs of G if every edge is turned on with probability p. These are 
polynomials in p. We have explicit recursive formulas in the case of K^, in which 
case the coefficients of / are integers. We can explore it numerically for others. 
The signature functions are certainly the same for isomorphic graphs but they do 
not characterize the graph yet. Two one dimensional graphs which are unions of 
disjoint cyclic graphs, we can have identical signature functions /, g if their vertex 
cardinalities agree. The union C5 U C5 and C4 U Cq for example have the same 
signature functions f^g. 

Since the global signature functions are not enough to characterize a graph, we 
can also look at the local dimension and Euler characteristic signature functions 

K (p) = Ep[l + dim{H, S{v)], K (p) = Ep [K{H, v)]], 

where H runs over all subgraphs and v over all vertices. Here K{H^ v) is the 
curvature of the vertex v in the graph H. Of course, by definition of dimension and 
by the Gauss-Bonnet- Chern theorem for graphs. 

hy{p) = f{p), ky{p) = g{p) . 

' ' vev vev 
Since for one-dimensional graphs without boundary, the curvature is zero every- 
where, also these do not form enough invariants for graphs. 

We could also look at the dimension and curvature correlation matrices 

A{p) = Corr[/i^(p), hy,{p)], B{p) = Corr[A:^(p), ku,{p)] 

which are nx n matrix- valued functions on [0, 1] if |V^| = n is the cardinality of the 
vertex set. We have 

^tr(A(p)) = Var^[dim], tr(5(p)) = Var^[x ] • 
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Figure 11. The dimension correlation matrix ^4(1/2) and the 
curvature correlation matrix B(l/2) for the Petersen graph G = 
P(3, 2) is shown. Below the dimension correlation matrix, we see 
the average local dimension. Below the curvature correlation ma- 
trix, we see the average curvature. The dimensions and correlations 
were computed in the probability space G{p) of all subgraphs of G 
which have the vertex set of G. The curvatures and dimensions at 
vertices which are not connected are uncorrelated. 




Figure 12. The dimension and Euler characteristic signature 
functions for the complete graph are 

f{p) = P{P^ -p^ -p^ ^ 5p^ - 3/ - 5p^ + lOp^ - 6p 4- 
g{p) = 5 - lOp + lOp^ - 5/ + p^^ . 
The bottom figure shows the case when the host graph is the cir- 
cular graph C7, where 

f{p)=p{2-p), gip) = 7{l-p) 

The function / is the same for every one dimensional graph, and 
in general g{p) = \V\{1 — p). 



4) 
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Figure 13. The left figure shows the logarithm of the num- 
ber of vertices having a given curvature on the probability space 
{V X G(7, 1/2)}, where V = {1,2, .., 7 } is the vertex set on the 
set G(7, 1/2) of graphs with 7 vertices. There are 47 different cur- 
vature values attained between —2 and 1. The right figure shows 
the pairs (dim(G), x(G)), where G runs over all 2(2) graphs of 7 
vertices. The Euler characteristic takes values between —5 and 7, 
the dimension of the graph takes values between and 6. 



Remarks. 



1. In the computation to Figure [13] we observed empirically that most nodes are 
fiat: 2784810 of the 7-2(2) vertices have zero curvature. Graphs are fiat for exam- 
ple if they are cyclic and one dimensional, which happens in 6! = 720 cases. The 
minimal curvature —2 is attained 7 times for star trees with central degree 6. 

2. Extreme cases of in the "dimension-Euler characteristic plane" are the complete 
graph with (dim, x) = (6,1), the discrete graph with (dim, x) = (0,7) as well as 
the case (dim, x) = (1,-5) with minimal Euler characteristic. An example is a 
graph with 12 edges, 7 vertices and no triangles obtained by taking an octahedron, 
remove all 4 vertices in the xy-plane then connect these 4 vertices with a newly 
added central point. 



1 2 


3 4 
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6 
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-19 -981 


-138043 


-4889125 


540429585637 


64 4096 


1048576 


1073741824 


4398046511104 



Table 3. The correlations of dim and x foi" until n = 7 for 
p = 1/2. These are exact rational numbers because dimension 
is a rational number and Euler characteristic is an integer. The 
correlation has become positive the first time at n = 7. For n 
larger than 7 we only numerically computed it with Monte Carlo 
runs in Figure [T4j 
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Figure 14. The left figure shows the expectation of the ran- 
dom variables dim and x G(n, 1/2) as well as Cov[dim, x] for 
n = 1, . . . , 14. The right figure shows the variance Var[dim], Var[x] 
as well as Gov [dim, x] on G(n, 1/2) again for n = 1, . . . , 14. While 
the variance of the Euler characteristic grows rather rapidly on 
that interval, the variance of dimension decreases on that interval. 
For n = 20, we measure it to be about 0.1. We can explore exper- 
imentally however only a very small part of the probability space 
G(20, 1/2) consisting of 2^2°) - 1.5 • 10^^ graphs. 
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min 
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1/2 





-1/2 -1 


-1.5 -2 


-2.5 -3 


-3.5 


-4 


-4.5 


max 


1 
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1 


1 1 


1 1 


> 1 > 1 


> 4/3 


> 3/2 


> 5/3 



Table 4. The maximal and minimal curvature which can occur 
at points of graphs with n vertices. This was settled forn = 1, . . . , 7 
by checking over all graphs. The minimal curvatures are obtained 
at star shaped trees, where curvature satisfies K{v) = 1 — deg('u)/2. 
For larger n, we ran Monte Carlo experiments over lO'OOO random 
graphs. 



Having looked at the random variables dim and x on G{n,p) it might be of 
interest to study the correlation 

Corr[dim, x] = E[dim • x] — E[dim] • E[x] 

between them. The extremal cases of size 0, order n graphs of Euler characteristic 
n and dimension or complete graphs with Euler characteristic 1 and dimension 
n — 1 suggest some anti correlation between dim and x; but there is no reason, why 
there should be any correlation trend between dimension and Euler characteristic 
in the limit n ^ oo. Like Euler characteristic, it could oscillate. 

While we have a feel for dimension of a large network as a measure of recursive 
" connectivity degree" , Euler characteristic does not have interpretations except in 
geometric situations with definite constant dimensions. For example, for a two 
dimensional network, it measures the number of components minus the number of 
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" holes" , components of the boundary dG the set of points for which the unit sphere 
is one dimensional but not a circle. While for geometric d dimensional graphs, 
x(G) has an interpretation in terms of Betti numbers, for general networks, both 
dimensions and curvatures varies from point to point and the meaning of Euler 
characteristic remains an enigma for complex networks. 
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